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Momentum and doping dependence of the static charge susceptibility x(q) i n the t — t' — J model 
is investigated. Correlations lead to a strongly momentum dependent renormalization of xCq)- 
The charge susceptibility near (ir, ir) region of the Brillouin zone is strongly suppressed as the hole 
density 8 is decreased. However, contrary to naive expectations, x(l) around q = (ir,0) and (0, ir) 
remains large and practically unchanged at 8 ~ 0.1 — 0.5. This effect is consistent with a tendency 
towards low-energy charge fluctuations with the wave vectors along the F — X direction, reported 
in earlier studies. Our main finding is that the above trends are amplified by J-driven pairing 
effects, indicating that the pseudogap formation may promote the charge inhomogeneity. The next- 
nearest hopping t' leads to weakening of the above momentum-selective renormalizations of x(q)- 
We analyze the effects of long-range Coulomb interaction, taking into account a layered structure 
of cuprates. As an application, the results are discussed in the context of bond-stretching phonon 
softening in hole-doped cuprates. In particular, a peculiar doping and momentum dependence of 
the electron-phonon coupling constant is found. 

PACS numbers: 74.72.-h, 71.27. +a, 71.10.Fd, 71.38.-k 



I. INTRODUCTION 

Low energy charge fluctuations and charge ordering 
becomes a hot topic in cuprates. Spatial modulation of 
the electronic states related to the local charge and/or 
bond ordering has been reported (see Ref. Q and ref- 
erences therein). Indirect evidence for the low-energy 
charge dynamics is obtained from phonon anomalies in- 
duced by hole doping in cuprates^^ These experiments 
motivate a theoretical study of the charge susceptibil- 
ity in correlated models. In general, one expects an 
overall suppression of the electronic density fluctuations, 
hence the related charge susceptibility, as one approaches 
the Mott insulating limit by removing the doped holes. 
On the other hand, it is also known that correlations 
may promote low-energy charge instabilities - e.g., so- 
called stripe physics in cuprates and other oxide mate- 
rials. These seemingly opposite trends indicate that the 
renormalization of the charge susceptibility by strong cor- 
relations is quite subtle process. 

Previous work on a charge response in t-J^^ 
and Hubbard models^ focused mostly on finite fre- 
quency charge response, x"(w,q), and on its frequency- 
integrated value, i.e. a structure factor N(q). These 
quantities provide an important information on electron- 
density fluctuation spectrum. In particular, Ref. pre- 
sented detailed calculations of x"(w,q) within a slave- 
boson framework. A nontrivial momentum structure 
of low-energy excitations has been found. Dressing of 
the doped-holes by underlying spin excitations - a phe- 
nomenon well known in the context of magnetically or- 
dered phase of t-J model - has also been captured within 
1/N expansion method for spin-disordered state. The re- 
sults of Ref. are in very good agreement with numerical 



Surprisingly, a static charge susceptibility Xq = x'i w = 
0, q) has escaped an attention. To our knowledge, no 
detailed discussion of the momentum and doping depen- 
dence of Xq m t-J model has thus far been reported. 
Meanwhile, this quantity which corresponds to the finite 
momentum compressibility contains an important infor- 
mation, e.g., about potential charge instabilities. The 
aim of this paper is to fill this gap. 

Specifically, we calculate x(l) m the t — t' — J model 
and discuss the doping and spin-pairing effects on x(q)- 
Consistent with known results, we observe that correla- 
tions may drive phase separation at small doping which is 
however eliminated by long-range Coulomb forces. Thus 
the charge instabilities, if any, are expected at interme- 
diate or large wave vectors. For large momenta, we find 
that the correlation effects are highly anisotropic in a 
momentum space. While x(q) is suppressed by a "hole- 
dilution" effect at certain parts of the Brillouin zone, it 
could even be enhanced for q along the T — X direction 
[from q = to (it, 0) or (Q,tt)}. As a result, a feature- 
less charge susceptibility x^Kq) ~ const of noninteract- 
ing electrons obtains a strong momentum dependence. 
Physically, these observations originate from a nontriv- 
ial momentum structure of low-energy charge excitations 
found in Refs.0IH3 

The main focus of the paper is to investigate how the 
above features in x(l) are affected by the J-term which 
induces a pseudogap in the fermionic dispersion. For- 
mally, this is done by considering fluctuations in the den- 
sity channel taking into account also the pairing fluctu- 
ations due to J-interaction. We find that the pairing 
effects cooperate with a Gutzwiller constraint and en- 
hance its momentum-selective renormalization of x(q). 
Another issue, raised in our study, is the influence of the 
next-to-nearest neighbor hopping t' , which is shown to 
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somewhat weaken the above anomalies in x(q). We also 
address a question how the momentum dependence of the 
compressibility is changed by the Coulomb interaction. 
We provide a realistic treatment of the Coulomb poten- 
tial accounting both for its long-range character and for 
the layered lattice structure of cuprates. 

There have been a number of discussions in litera- 
ture on how the correlations renormalize electron-phonon 
coupling, 11 ! 12 ! 13 ! 14 ! 15 . 16 ! 17 ! 18 In Holstein-Hubbard type 
models (relevant to the problem of oxygen vibrations 
coupled to the electron-density) it was found that a pe- 
culiar "forward scattering feature" may develop due to 
correlations. We will discuss a connection between this 
observation and our findings for the charge susceptibility 
x(q). Related to this issue is the bond-stretching phonon 
anomalies in cuprates, which are discussed in the last part 
of this paper. This part extends the previous studyi^^ 
of the phonon-softening problem by including the pairing 
and t' effects. 

The rest of the paper is organized as follows. Section II 
describes the formalism and discusses x(q) in the i-only 
model. Sections III and IV focus on the effects of J and 
t' terms, correspondingly. In Section V, we derive a mo- 
mentum dependence of Coulomb potential in the layered 
lattice structure, and calculate x(q) at presence of these 
interactions. The last Section VI discusses renormaliza- 
tion of the electron-phonon coupling by correlations. 



II. THE MODEL AND FORMALISM 

The t-J Hamiltonian is 

ij i 
+ J ^(SiSj - -UiTlj) + ^ VijUiUj. (1) 

One of the useful approaches in treating the local con- 
straint on the fermion occupation number, ^ cj^c^ < 
1, is the slave-boson representation Cj CT = fi a b\. The 
above inequality is then replaced by the constraint b\bi + 

E CT ft, fi° = N / 2 > with the Physical case of N = 2. It 
is convenient also to consider the limit of the large num- 
ber of spin indices (flavors), N 1, since calculations 
are simplified this way (see, e.g., Ref. and references 
therein). 

Formally, the slave-boson approach implements the 
constraint of no double occupancy at the operator level. 
In the large- N limit, the mean number of bosons (bjbi) is 
large and bosonic amplitude contains a large c— number 
component, a bosonic condensate. It was noted however, 
that the phase of each slave boson is a gauge degree of 
freedom, which is eliminated from the action by promot- 
ing the local Lagrange multipliers into time-dependent 
fields. The remaining degree of freedom for slave-bosons 
is their real- valued amplitude, bj 



can formulate the "radial-gauge" representation for slave- 
bosons. 

Below we show how the leading-order results obtained 
for the charge susceptibility in the large-iV slave-boson 
approach 5 * 7 " can be reproduced in a simple manner. To 
clarify our approach, we let J = first, so that only a 
hopping term is present. We represent the constrained 
fermions as 



(2) 



with rii the fermion occupation number, rii = c ia Ci„. 

We are interested in the charge susceptibilty of the 
system, therefore we consider small fluctuations of the 
fermionic density, Uj, around its uniform equilibrium 
value, n. We expand the hopping term up to the sec- 
ond order in 

</>,: = Sni/2 = [rii — n)/2 , 
with the factor 2 introduced for later convenience: 
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2(1 -n) 



(3) 



In the Fourier representation the first term above be- 
comes a fermionic dispersion with the renormalized am- 
plitude, the second term describes the scattering of the 
fermions on the fluctuations of density. The third term 
has more complicated structure, for our purposes it is 
enough to consider its part, which is diagonal in fermionic 
momenta, cx cL a c^ a with k\ = k 2 . Then we write 



Hi 



Cfco 



(tk + t k +q)c{ a C k+q , a (j)q 

k.q,a 



kqa 



— c' ka C ka (j) q (j)- q . (4) 



with 



& = -(l-n)tk+H, (5) 
tk = 2£(cos k x + cos k y ) — At' cos k x cos k y , (6) 

henceforth we set t = 1. The constraint for n, to be a 
number of on-site electrons leads to the appearance of 
the local Lagrange multipliers in the action 

-fl q [2(j)^ q - J2 k a C ka C k+q,*) , 

so that the scattering term takes the form 



Tj and one 



k,q 



(7) 
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with Uk q = 4> q {tk + tk+q) + AV The average number of 
fcrmions is defined by n — 2^ fc 7ij?(£fc) with np(x) the 
Fermi factor. One can integrate out the fermions now, 
and obtain the effective low-energy action: 



Ft - E [^"9 ( G 



W 



Or. 



20,^_,(ni + 1) - /i g M- g n ] , (8) 



(9) 



E % ± %^ (^&)-Ma +g ))- (10) 



Requiring the zero variation, bT t j§\iq — 0, we determine 
the values of Lagrange multipliers 



nq = -0 9 (i + ni)/n o , 



(ii) 



Inserting these values into (JBJ, and recalling that <f> q — 
riq/2 at q ^ 0, we find 



E —r~ K ~ n 



(i + ni 
n n 



The above equation for the free energy in the harmonic 
approximation should be compared to the general expres- 



sion T = J2 Q ( 2 Xq) l n q 



with the static charge 



susceptibility Xq- The value of \ q as determined from 
(|12[) coincides with the result of Ref. in the limit of 
uj = [however, a so-called polaron correction due to the 
higher 1/AT-termi is absent in Ea. lfT^l) ]. 

Let us first qualitatively analyze the above formula for 
the charge susceptibility. We introduce the doping level, 
or the concentration of holes, 

5 = 1 — n, 

and extract this factor from the dispersion, £ q , and the 
chemical potential, /z, 



(13) 



Note that jl is positive and proportional to the doping 
level at t' = 0. 

After some straightforward rearrangements of Eq. (|12|l , 
the expression for x(l) can be represented in the follow- 
ing form: 



Xq = x^ 



(6 - a q ) 2 + 2{fi5 + /3 q )^ 0) 



(14) 



Bare quantities Xq°^ > c^q and /3 q are given in terms of the 
" noninteracting" dispersion, £jj, ^ as follows : 



i 



AO) _ AO) 



Mt^-Me&J), (is) 



(0) 



corresponding to the bare susceptibility, and 



AO) 

^k+q 

k 



AO) _Ao) 



(0), 



E 



t(0) t(0) 



(^(^ 0) )-n F (^ g )).(17) 



MO' 



(0) 



(16) 



t(0) 



(0) 



One can show that the function /3 q in i|17fl is posi- 
tive. Eci. (|14fl should be understood as a renormaliza- 
tion of the charge susceptibility by correlation effects, 



Xq/X 



(0) 

q 



Gq(S), where a function G q (<5) is given by 



a fraction in Eq. l|14|) . This momentum and doping de- 
pendent factor results from the action of the Gutzwiller 
constraint in the density channel. Whereas a noninter- 
acting susceptibility Xq^ is a featureless function in the 
absence of nesting (at 8 = 0, t' = 0), the correlation ef- 
fects bring a pronounced momentum structure in Xq via 



the function G q (<5). 

At large doping levels, G q ((5) eventually approaches 
unity. The action of this factor at small doping is highly 
momentum-selective. Inspecting Eas. 11141171 at t' = 0, 
(12) one nn ds that G q (<5) oc 6 for q ~ (tt,tt) (with omitted 
logarithmic corrections). This is simply understood as a 
reduction of density fluctuations due to a removal of the 
holes. (Alternatively, one may say, that the checkerboard 
structure in positions of small amount of doped holes is 
the least energetically favorable.) At small momenta, 
however, the effect is opposite and one has G q ((5) oc 1/5. 
This means a divergent compressibility as one approaches 
the Mott limit, in accordance with previous studies of 
HubbardSi and t-J modelSj 22 i 23 ' 24 and reflects a well- 
known tendency towards phase separationiSSiSSiSi Com- 
petition between these two effects - a hole dilution and 
phase separation - leads to a nontrivial momentum struc- 
ture. It is interesting to note that this structure is com- 
plementary to that in the spin sector, where correlations 
enhance the spin susceptibility at q ~ (it, ir)^^ but not 
at small momenta. 

These qualitative observations are further illustrated 
by the numerical calculations. In Fig.^ we show the mo- 
mentum dependence of x q along the symmetry lines in 
the Brillouin zone. The above behavior of Xq with doping 
is clearly visible at the symmetry points. In order to em- 
phasize this nontrivial momentum dependence induced 



by the Gutzwiller constraint, we plot Xq/X 



10) 



Gq(^) 



in Fig. H for several dopings 6 ~ (0.1,0.2,0.3,0.5). To 
see the doping dependence in more detail, we show x q in 
Fig. |3 as a function of doping at three symmetry points 
q = 0, (7r, 0) and (tt,tt). One observes that the curves 
for x q ^ eventually turn down at small doping 6. Re- 
markably, the value of Xq at (tTj 0) upon decreasing S 
is somewhat enhanced before the downward turn, which 
shows the competition between two trends: phase sepa- 
ration and hole-dilution. 

Small momentum anomalies are eliminated in reality 
by Coulomb repulsion. A detailed study of this problem 
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FIG. 1: The inverse non-uniform compressibility calculated 
for tight-binding spectrum. 
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FIG. 2: Renormalization of the susceptibility due to the 
Gutzwiller constraint. 



is presented in Section V. We show there that a pro- 
nounced momentum structure of x(q) with very differ- 
ent doping dependence at (it, it) and (ir, 0) regions still 
remains at the presence of Coulomb interactions. 



III. J-TERM: PSEUDOGAP EFFECTS 



We consider now how the above observations change 
at the presence of pseudogap effects induced by J-term 
in the Hamiltonian. In the spirit of large- N slave-boson 
theories, we refer to pseudogap as a fermionic gap arising 
from mean-field decoupling of the superexchangc interac- 
tion in the pairing channel. A question addressed here is 
that how such a gap and fluctuations of the pairing field 
around uniform mean-field solution will affect the charge 
susceptibility. 

The four-fcrmion J-term can be represented in the 




FIG. 3: Doping dependence of the charge susceptibility (up- 
per panel) and its inverse value (lower panel) at the symmetry 
points. 



form 



n 



J = —~ ^2 C !l,T C fe2,l C fc3,iCfc4,T(^fc4-fcl + Jk3-kl), 



fel,fc2,fe3 



(18) 

with fej = fci + &2 — ^3 an d the nearest-neighbor interac- 
tion Jk = 2J(cosk x + cos k y ). Introducing the quantity 



with 



7 fc = (cosk x ±cosk y )/2, 



we represent the J— term as 



nj = -47 J2 

q.a— ± 



(19) 



(20) 



(21) 



This expression can be decoupled by the Hubbard- 
Stratonovich transformation as follows 



q,± 



Id* I' 



(22) 



Here d~ and <i + stand for the amplitudes of the d— wave 
and extended s— wave pairing, respectively, in the chan- 
nel with non-zero total momentum, q. 
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The total Hamiltonian is then quadratic in fermions, 
which interact with the fluctuations <fi q ,n q and d~ . 

We assume that the d-wave pairing sets in, which cor- 
responds to the non-zero value of the order parameter 
o?q = A. The spectrum is given by er| = ^ + A?, where 
Afc = dyj^. The self-consistency gap equation in the 
small-coupling limit reads 



1 = 2J yi2£l! tanh £^ 



(23) 



Restricting our consideration by the quadratic terms 
in the above non-uniform fluctuations, we would like to 
obtain an expression similar to 1|12|1 . but in the presence 
of the pairing. After some standard analysis, we arrive 
at the bosonic-type action of the form 



9 

{<f> q ,H q ,d+,d-), 



*9 



(24) 
(25) 



with the matrix 



M q = 



/ uj q -U 2 


-i-n x 


At 


AI \ 


-1-Hi 


-Ho 


At 


A 


At 


At 


{U)~ x -D ss 


-D sd 


V 


^0 


-D sd 





Here in the low-temperature limit 



(26) 



tk — tk+q £>. - Li, 

k 



n Ek 

e+e- - e+C- + A + A_ 



X>++*-) r 



2e + e_(e + + s- 



A + g_+g + A_ / 7 + 
2e + e_(£ + + e_) ^ 



(27) 



(28) 



and 




(29) 



and we used the shorthand notation £± — £k±q/2i £ ± = 
£k± q /2, etc. 

We are interested in the charge susceptibility, which 
should basically be determined by integrating out the 
pairing fluctuations in (|25|) and setting the Lagrange mul- 
tipliers /i q to their saddle-point values. Both steps are 
essentially the same for the quadratic action, so that the 
needed compressibility is given by the upper left element 
of the inverse matrix M q , namely 



X, = 2[M- 1 ] U . 



(30) 



Let us briefly discuss here how the inter-site Coulomb 
interaction Vn is included in our formalism. We write it 



in the form Vij ( n i ~ n)(n 



2 J2 q Vq't'q'f'-qi an d see that this interaction modifies the 
only matrix element in Eq. (|26|l . so that 



[M, 



q\ll 



^ + 2Vg 



(31) 



The last equation shows that we treat the interaction V q 
within the RPA scheme. It can also be shown that 



X V =0,q + V 9> 



(32) 



i.e. one can calculate Xq at V q = 0, and include V q ^ 
afterwards. At the same time, the inclusion of V q ^ 
should be done from the beginning for the calculation of 
d— wave susceptibility, Xdd(q)> see below. 

It is worth also noting that in the absence of pairing, 
Afe = 0, the coefficients D a p in 126f) remain finite, and 
vanish. It means that in the harmonic approximation 
the superconducting-type fluctuations affect the charge 
susceptibility only at finite A. At the same time, in more 
general treatment, the superconducting fluctuations d 
affect the density fluctuations in the higher orders even 
at Afe = 0. Considering multi-tail fermionic loops, one 
obtains, e.g., the terms in the bosonic action of the form : 



a; 



'cid k2 d k3 



etc., where the coefficients A\ , B are defined by the 
Feynman diagrams shown schematically in Fig.^Ja) and 
Fig.0Jb,c), respectively. One can notice that the diagram 
in Fig.QJb) corresponds to an analog of Maki- Thompson 
contribution to paraconductivity and the diagram in Fig. 
[Uc) reflects the "density of states correction", see Ref. 
l29t The contribution of the diagram Fig. Ufa) is zero at 
Afe = 0. However, in the presence of the pairing conden- 
sate, one external bosonic field, d~ , sets to a constant, 
contributing in the lowest order to the term A^ in 12611 . 
etc. The calculation of such fluctuation corrections to 
the charge susceptibility is clearly beyond the scope of 
the present study. 

In our numerical calculations we set J = 0.3, and de- 
termined A from the self-consistency equation l|23l) at 
T = 0. One should note, that at T = Eq. (J23) has a 
solution A ^ for any doping, although A can be ex- 
ponentially small. In our calculation, we ruled out the 
solutions with |A| < 10~ 3 i, thus implicitly setting the 
temperature to be very small but finite, T ~ 10~ 3 t. The 
matrix M is then found from J2ZJ, (PSJ, and \ q 

from (|2Dl. 

The obtained results are shown in Fig. |SJ Compar- 
ing the upper panel in Fig. [3] with the previous Fig. 
n w e observe that \ q is still rather flat in the whole 
Brillouin zone at large dopings. At smaller dopings the 
above momentum-selective features are enhanced by the 
pseudogap appearance, with somewhat increase of x^ 1 
at q = (jr, 7r) and a new qualitative change at small 
wave-vectors. Namely, the compressibility attains neg- 
ative values at finite S ~ 0.15, which means that the 




FIG. 4: Multi-tail fermionic Feynman diagrams, leading to 
higher-order terms in the bosonic action. Fermionic Green 
functions shown by lines with arrows, bosonic fields cj>, d~~ 
stand in the vertices. For simplicity we do not show the in- 
ternal momenta. 



Gaussian action ~ Xq(^ n q) is unstable in rather ex- 
tended range of dopings and the analysis of the next 
orders in 8n q is needed. This instability is accompa- 
nied by the divergence in the d-wave susceptibility, Xdd, 
defined as the ^[M" 1 ^, as is seen in the lower panel 
in Fig. [SJ We remind that the latter quantity is al- 
ways positive Xdd ~ ^ 2 'A m the absence of feedback 
from density fluctuations to the superconducting ones, 
i.e. when = 0. The unstable pairing part of the ac- 
tion ~ ^ Xdd (Q^dq \ 2 particularly means that the uni- 
formly paired ground state determined by the gap equa- 
tion JZ| is no longer justified. 

Our finding that J-pairing fluctuations and charge fluc- 
tuations grow up concomitantly could be understood as 
a dymanical modulations of pairing amplitude consistent 
with the results by Vojta et al£&Q It is also noticed that 
a dramatic enhancement of the charge susceptibility at 
small momenta due to J term is consistent with previous 
reports (see, e.g. Fig. 3 of Ref l2rj) that the superex- 
change interaction increases a tendency towards phase 
separation. Close to such instabilities, the higher order 
terms (beyond the Gaussian action) should be included 
in the theory, which problem deserves a separate study. 



IV. NEXT-NEAREST HOPPING 

We discuss now the effects of next-to-nearest neigh- 
bor hopping t' which is present in cuprates and has in 
fact been suggested to be a key empirical parameter 
for superconductivity^ Effects of t' on physical quan- 
tities such as spin and fermionic excitations in the t-J 
model has been found to be substantial, see, e.g. Ref. 
l33l Concerning the charge compressibility, several stud- 
ies found that t' hopping reduces a tendency towards 
phase separationiSSiMiS 

Consider first the qualitative effect of t' 7^ in the 
absence of J term. Using l|14jl. one can still show that 
the renormalization factor G q (<5) ~ S at q = (tt,tt). At 
the same time, for t' ^ the chemical potential jj, does 
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FIG. 5: Behavior of inverse charge (upper panel) and d-wave 
pair (lower panel) susceptibilities at J = 0.3. 



not vanish when 8 — > and we have finite Gq(5) ~ at 
small momenta. It means that finite (positive) t' reduces 
the tendency to the phase separation at small doping, 
consistent with previous work. It is interesting to note 
that the negative t' has an opposite effect. Indeed, in 
this case the chemical potential /Et is also negative. By 
inspecting Eq.jnj one observes that this may lead to a 
negative values of the susceptibility at small momenta 
and doping, indicating an instability of the uniform state 
for t' < case at small dopings. 

We showed in the previous section that the inclusion of 
J term drives the system closer to the instability point 
for charge fluctuations at t' = 0, and the same thing 
should happen when next-to-nearest neighbor hopping is 
present. To verify it, we recalculated x q for t' = 0.3, J = 
0.3 for the same values of doping 5 as above. The results 
are shown in Fig. El Comparing it to Fig.[Sl one confirms 
that the finite next-nearest hopping t' > Q somewhat 
stabilizes the charge fluctuations. The comparison to Fig. 
n shows, however, that the effect of the J term is still 
dominant and \q is (nearly) divergent at 5 = 0.1. 

We should emphasize that the above statements on the 
role of t' and J terms do not explicitly rely on the one- 
particle properties of the spectrum, such as van Hove sin- 
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FIG. 6: Inverse charge susceptibility at t' = 0.3 and J = 0.3. 
Anomalies at small momenta seen in Fig. 5 are suppressed. 



Here, e a b and e c are zero-frequency dielectric constants 
determined usually from optical data, and q 2 = q 2 + q 2 . 
(a and b directions are assumed to have the same e). We 
use a notation Q = (q, q z ), with in-plane and out-of- 
plane components, q and q z , respectively. For q z = 0, 
the above equation gives V(q) = Aire 2 /e a bq 2 , a conven- 
tional 3D potential. We recall that in 2D (the case of 
infinitely separated planes) V(q) oc 1/q, and we discuss 
the crossover between these two regimes below. 

We argue here that for our analysis it is possible to 
neglect a momentum dependence of the dielectric con- 
stants, because at low energy they are mostly contributed 
by (dispersionless) optical phonons and nearly localized, 
high-energy electronic processes. In this case the real 
space representation of (|33|) . valid up to interatomic dis- 
tances, reads as follows: 



gularities and flat parts of dispersion around (0, n) points. 
Our discussion includes two-particle Green's functions, 
both particle-hole and particle-particle fermionic loops, 
and the eventual integration over bosons in the effec- 
tive action, i.e. obtaining Eq. (|3U|I from Eq. I|24|) . corre- 
sponds to simultaneous resumming of RPA series in both 
Gutzwiller and J-term channels. 

Summarizing here, the static susceptibility shows no 
structure in the Brillouin zone at large value of the hole 
doping, S ~ 0.5. This could be expected for a system 
without strong correlations and with a large Fermi sur- 
face, since in the 2D Fermi gas \q — const at q < 2kp. 
The flat shape of Xq at large S is rather insensitive to the 
values of the second hopping and pairing. At smaller dop- 
ings, Xq demonstrates a pronounced structure in (/-space; 
the tendency to long-scale phase separation is somewhat 
weakened by finite values of the second hopping, but the 
pairing fluctuations, induced by the J term, dominate 
and eventually make the system unstable both in charge 
and pairing channels. 



V. LONG-RANGE COULOMB INTERACTION 

Charge susceptibility is strongly influenced by a non- 
local repulsion between the holes. Quite often (in nu- 
merical studies in particular) these interactions are ap- 
proximated by a nearest-neighbor potential Vi, which is 
already sufficient to observe the suppression of phase sep- 
aration effects discussed above. We are however inter- 
ested in a more detailed momentum dependence of Xq- 
For this purpose, one has to use more realistic, i.e. long- 
range form of the Coulomb potential. We consider first 
its momentum dependence in a layered cuprate structure, 
taking into account a discrete nature of the lattice within 
the planes as well. 

At small momenta, where no lattice structure is rele- 
vant, a continuum limit applies: 



Vb(Q) 



Aire 2 



e a bq 2 + e c q£ 



(33) 



V C (R) 



1 



y^aWc [(e ab /e c )z 2 + r 2 ]V2 



(34) 



Here, r is a distance within a6-plane, and R 2 = r 2 + z 2 . 
In isotropic case, e a b = e c , a familiar expression e /eR 
follows from this equation. 

Let us consider now a lattice with periodicity a within 
the planes and d along the c-axis (in La2Cu04-structure d 
is a half of the c axis lattice parameter, i.e. d = c/2). We 
determine the Coulomb repulsion Vij between the elec- 
trons, referring to the sites i and j. In our tight-binding 
situation the electronic wave-functions are almost local- 
ized around the i-th ion and their amplitude squared 
gives the density around this ion. We denote this density, 
or charge distribution function, by /(R — Rj), and write 

Vij = J dTL'dR." f(TL' - R)/(R" - R;)V C (R' - R") 
rfQ 



-\f(Q)\ Z V c (Q)e 



■Q(Ri-Rj) 



(35) 



where the integration is over the whole continuum and 
Vc is given by Eqs. I|33|l and l|34|) in a momentum and 
real spaces, respectively. The Fourier transform of Vij 
then reads as 



V(Q) = J2 V ii eiQ(Ri ~ ni) 



(36) 



l - |/(Q + G 3 )| 2 Vb(Q + G 3 ) - Vu 



a 2 d 



with G3 the 3D wave vector of reciprocal lattice. 

On the physical grounds, one expects that a doped 
hole (the Zhang-Rice singlet) is a rather extended object 
in the afr-plane and nearly localized in this plane. A 
reasonable choice for a hole-shape function is thus 

/(R) = (K 2 /2n)e- Kr S(z) 

which decays at distances 1/k in the plane. Physically, 
the size of the Zhang- Rice singlet should at least be about 



Cu-0 distance, so k ~ 2/ a might be a representative 
value. A momentum counterpart of the latter function, 



/(Q) = (i+? 2 /«r 3/2 



(37) 



should then be understood as a formfactor of the Zhang- 
Rice singlet. 

Given that the formfactor /(Q) is independent of the 
q z component, the summation over G z = 2im jd (n — 
0,±l,±2,...)in l|37[) is easily performed for any q z . The 
result for our primary case of interest, q z = 0, is 

V q = V(q, 0) = + G 2 ) |V (0) (q + G 2 ) - V H (38) 

G 2 

where G 2 = 27r(n, m)/a is the reciprocal wave vector for 
square lattice and 



y(°)(q) = 

V = 



v 



qa tanh((7/qo) 
27re 2 



q = 2/d={2/d)y/e c /e ab 



(39) 

(40) 
(41) 



Here d is an effective interlayer distance. The poten- 
tial (|3l^|l interpolates between 3D Vqo/q 2 a and 2D V/qa 
limits at small q <C qo and large q ^ qo momenta, re- 
spectively. This crossover at go reflects the fact that the 
planes are independent at large momenta. For La 2 Cu04 
compound where e c /e a b ~ 1/2)2* the value of crossover 
momentum is estimated as qo — OS/a. 

The function l|38[) is explicitly periodic in q— space. In 
the particular model for the formfactor, Eq. I|37[) . the 
subtracted term Vu is evaluated as Vu — 3naV/32 . In 
general, V q is sign-reversal function in the Brillouin zone, 
as it should be in view of £ Q V(Q) = V(R = 0) = 0. 
In case of a continuum limit for the planes^ (instead 
of the tight-binding model used here), the Coulomb po- 
tential would be given solely by Ea. (|39|l . the result for- 
mally obtained from Eq.JSEl by setting Vu = 0, /(q) = 1 
and taking G2 = term alone. Finally, we note for the 
completeness that V(q, q z ) for arbitrary q z is again given 
by Eq.® but the potential V^ 0) (q) for a contimuum 
limit in Ea. H39|) must be replaced by V' '(q, g 2 ) = 
VW(q)/(l+F%), where F z = sin(g z d/2)/ sinh(gj/2)^i 

At moderate k ~ 1 and small momenta, q < k, (or, 
alternatively, at distances exceeding the size of Zhang- 
Rice singlet) only G 2 = term in (|38|l contributes and 
V q ~ V"(°)(q). The rapid decay of |/(q)| 2 , E3, cuts off 
the values of V q at larger momenta, q > n. 

The opposite limit, > 1, corresponds to the point- 
charge approximation and the sum in (|38|l formally di- 
verges with k. This divergence is cancelled by the above 
on-site term Vu oc k. In Appendix, we provide another 
representation for V q in the point-charge limit; it works 
increasingly well for na > 3. 

Fig. shows q-dependence of the Coulomb potential, 
Eq. (l38|l . as function of parameter k at fixed qo — 0.8 /a. 



> 

C 

> 



K= 1.0 

k = 2.0 
k = 4.0 
short-range V 




FIG. 7: Coulomb potential Vq/V for different values of k. A 
crossover momentum go = 0.8/a. For comparison, a dashed- 
dotted line shows a simple model which includes NN- and 
NNN-repulsion only. Dotted line is calculated from Eq. l39H 
neglecting the lattice structure within the planes. 



Vq/V evolves from the point-charge limit with visible 
sign-reversal character to nearly positive curves as k de- 
creases. For comparison, we also show a frequently used 
simple short-range repulsion model, 



V q = 2^(008^ 



cos q y ) + 4V 2 cos q x cos q y , 



with Vi = V/2tt, V 2 = V1/V2 (dashed-dotted line), and 
the result that would be obtained when lattice structure 
is discarded within the planes (dotted line). The latter 
is always positive as should be for the charges in a con- 
tinuum. 

Let us turn now to the charge susceptibility and con- 
sider how it is influenced by long-range Coulomb interac- 
tions. First, we estimate the energy scale V in Ea. (|40(l . 
Using a representative value ^/e a (,e c ~ one finds 

V ~ 0.8 eV, which is about 2 in units of t~ Second, 
we should in principle complement the Coulomb poten- 
tial with a short-range interactions between the holes, 
stemming from local physics. One such a contribution 
is that of well-known " missing J-link" , which gives NN- 
attraction of the scale of J(siSj — 1/4) ~ — (0.1 -j- Q.2)t. 
Yet another local interaction is mediated by the bond- 
stretching vibration of an oxygen shared by the two NN- 
holes. This contribution is repulsive at low energy limit 
because of the coupling geometry (see for details the next 
section), and is given by a half of the polaron binding en- 
ergy Eb/2 ~ t/A (estimated below from the phonon shift 
induced by doping). Altogether, these two local NN- 
contributions of different sign tend to cancel each other 
and their small net result could be neglected. Hence, we 
focus on the Coulomb repulsion. 

Fig. [S] shows a momentum dependence of the charge 
susceptibility at the presence of Coulomb interactions. 
The parameters used: qo — 0.8/a and V = 2. Com- 
pared with a pure t-J model result in Fig. 4, one observes 
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FIG. 8: Behavior of inverse charge susceptibility at J = 0.3, 
t' = 0. Coulomb repulsion with V — 2, k — 2/a and go = 
0.8/a is included. Compare with figure 5. 



that Vq eliminates phase separation effects. We find also 
that small momentum divergences of the pairing fluctu- 
ations are suppressed, too. However, Coulomb interac- 
tion effects are not significant at larger momenta. All 
the local correlation effects, which lead to a pronounced 
anisotropy of the charge susceptibility and its nontrivial 
doping dependence along the T — X direction, remain 
intact. As expected, the main effect of Coulomb interac- 
tions is to move a potential charge instabilities to a finite 
momenta, as seen in Fig. [H] The divergence of Xq 1 a t 
q — > corresponds to x(q = 0) = 0, a well-known screen- 
ing phenomenon in the presence of long-range Coulomb 
repulsion. 



VI. APPLICATION: PHONON SOFTENING IN 
CUPRATES 

Density fluctuations determine doping induced phonon 
renormalization, as discussed earlier in a framework of 
slave-boson methodji^iS . and also by an exact diagonal- 
ization of small clusters^ In particular, a broad and 
anomalous lineshape of the bond-streching phonons with 
momentum at (tt, 0) direction has been found^S while 
no such anomalies were present for (tt, tt). The effect 
is strongly doping dependent^ We consider now how 
the phonon renormalization effects are modified when the 
pairing (pseudogap formation) is included. Specifically, 
we address here a phonon softening problem and discuss 
the results in the context of experimental reportsi^*^ 

A doped hole couples to the bond-stretching vibrations 
of its four oxygen neighbours^: 



H = g 



E 



m(u' 



In a momentum space this reads as 

H = i\J E b ujQ sin(q a a/2)(a a _ q + a\ c 



(42) 



(43) 



Here, a = x and y denotes polarization of the oxygen 
displacement, and u>q = y/K/m is the phonon frequency 
determined by spring constant K and the oxygen mass 
m. Electron-phonon coupling strength is conveniently 
quantified by a binding energy 



E b = 2g 2 /K, 



(44) 



that would be gained in case of a static hole. By fitting a 
slave-boson theory to the experimental data on phonon 
softening and linewidth in cuprates, an estimation Eb ~ 
t/2 was obtained in Ref. lift 

We assume that the energies of density fluctuations are 
higher than the phonon energy — it should be valid not 
too close to charge instability. Within this assumption, 
we can use our static Xq to estimate phonon softening, 
Wq(<5), which is obtained as follows: 

£*Xq(l-7q). (45) 



8lo d 



LOO 



^ =l-y/l -2£ 6Xq (l 
LUq * 



7q) 



where 7 q = (cosq x + cosq y )/2 and Xq includes the 
Coulomb repulsion. Note also, that without correlations 
the product -E&Xq — Eh/At. One can therefore introduce 
a dimensionless quantity, A'°) = Eb/4t, which could be 
regarded as a "bare" coupling constant in the problem. 
According to Ref. [Hj, A' ) ~ 1/8, justifying a perturba- 
tive treatment. Eq. (|45|l reads now as 



(46) 



Swjua = A q = A(°)4i Xq (l- 7q ). 



a . q 



Factor 7 q in Eq. l|4()|) . which stems from the coupling 
geometry, would suggest a strongest softening for the 
full-breathing mode, that is at q = (tt, tt). However, 
strong correlations change a momentum dependence of 
Xq dramatically, by suppressing it at (tt, tt) and enhanc- 
ing around (tt, 0) points. One may say that correlations 
lead to the redistribution of the effective electron-phonon 
coupling in a momentum space. As a result, softening 
becomes strongest at (n, 0), consistent with experiment. 
This explanation of Ref. [19| is further supported by the 
present calculation, including the pseudogap and t 1 ef- 
fects. 

Fig.|5|shows a general form of the renormalization fac- 
tor A q /A (0) along particular directions in the Brillouin 
zone. This figure is in obvious correspondence with the 
above findings, and particularly, with Fig. 

Fig. ^5] presents more detailed doping dependence of 
Atj/A^ ) at q = (tt, tt) and q = (tt, 0). When multiplied 
by a bare constant X^°\ these curves correspond to the 
phonon softening Sujq/uJo. As the latter is about 15 - 
20% for (tt, 0) phonon in optimally doped cuprates;^ a 
bare constant of the order of 0.15 - 0.20 is required 
to fit the observed data^i. A striking similarity with the 
observed doping dependence^ is worth to be pointed out 
here: both in experiment and in our theory phonon soft- 
ening along (tt, 0) direction is almost independent on dop- 
ing in a wide region above 8 ~ 0.12. While such a trend 
was already found earlier j2S J-pseudogap effects dramat- 
ically enhance the charge susceptibility along (tt, 0) at 
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FIG. 9: Momentum dependence of the effective coupling con- 
stant (in units of the bare coupling A' ') which determines 
renormalization of the bond-stretching phonons. Parameters 
used: t' = 0.3, J = 0.3, V = 2 and k = 2/a. 
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FIG. 10: Renormalization of phonons (in units of the bare 
coupling constant A' ') at the symmetry points as a function 
of doping. The parameters used are as in Fig. [§] 



small <5, hence it obtains nearly flat doping dependence 
rather unexpected in view of hole-dilution physics. 

For further comparison of our theory with the available 
experimental data, we plot a momentum dependence of 
Aq/A^ ) along r — X direction for several values of doping 
in Fig. ^2 One finds that visible deviations from a sim- 
ple cosine curve increase at smaller dopings, in general 
agreement with experiments^ 

While present calculations do capture the most anoma- 
lous experimental findings - stronger and nonlinear dop- 
ing effects along (tt, 0) direction - a quantitative com- 
parison is much less satisfactory. In particular, rather 
sharp kink-like change in doping dependence is observed 
in experiment at about S ~ 0.12^2^ while it is found in 
our theory at lower doping. One obvious reason for this 
discrepancy is that, strickly speaking, we cannot quan- 
titatively address the phonon softening problem by us- 
ing our static charge susceptibility. This is because the 
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FIG. 11: Renormalization of phonons along the (0,7r) di- 
rection (in units of A' - 1 ). The parameters used: t' = 0.3, 
J = 0.3, V — 2, k = 2/a and qo — 0.8 /a. A simple cosine 
curve is shown for comparison by dotted line. 



charge fluctuations for momenta along (7r,0) direction 
extend to low energies (comparable to phonon ones), as 
shown both in a slave-boson theorjii and in the numerical 
work^iS A dynamical susceptibility is therefore required 
which is however beyond the scope of present work. Yet 
another reason is that, focusing mainly on the pairing ef- 
fects, we did not include an effective hopping contribution 
stemmimg from "Fock" decoupling of the J interaction. 
This contribution, which is a fraction of J/ 1 jil will sta- 
bilize a Gutzwiller band narrowing effect at doping levels 
~ (0.2 — 0.3) J jt below which a linear doping dependence 
Xq oc 5 due to a hole-dilution effect sets in. This is ex- 
pected to shift a kink feature to higher dopings as in 
experiment. Further, a quantitative description would 
also require the inclusion of spin-polaron effects beyond 
the leading 1/N approximation 7 ^, and a coupling of the 
charge fluctuations to a collective spin modei 4 ^ 

In general, it seems that a kink feature in the doping 
dependence of (tt, 0) phonon softening^, which appar- 
ently parallels with the so-called "Yamada plot" for the 
spin incommensurability 43 , provides an interesting test 
case for theory. Its initial linear behavior oc S is expected 
and can be explained in terms of sum rule arguments. 17 
The saturation above certain doping level is well cap- 
tured by slave-boson theories, but it is not fully clear at 
present why such an abrupt regime change happens at 
around the "magic" doping ~ 1/8 concomitant with the 
saturation of spin incommensurability. 

Finally, it is interesting to notice that the renormaliza- 
tion of density fluctuations and the so-called charge ver- 
teWfc, q) of Refs lIlTal (denoted by T(p, q) in RefsE 
mil have a common origin. Reflecting this, the charge 
susceptibility can readily be expessed via j(k,q) (see, 
e.g. Eq.(16) of Ref If3T) . In fact, in case of zero-frequency 

and small momenta one finds 7(fc/,q) ~ 5 ■ (Xq/Xq^)- 
This relation makes it clear the origin of strong momen- 
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turn structure in electron-phonon vertex function found 
in Refs lT2TITM3l5lll6Jll8l — this simply reflects highly 
momentum-selective action of the Gutzwiller constraint 
on density fluctuations as we emphasized in previous sec- 
tions of this paper. Indeed, as Xq/Xq K l/<5 at small 
momenta, one realizes that the effective electron-phonon 
interaction is essentially the bare one. At large momenta, 
however, it is strongly suppressed as the hole density S 
is reduced^ This results in a "predominantly forward 
scattering" of electrons on phonons in cases when this 
coupling is located in a density channel (which is the 
case for the scattering on bond-stretching phonons). In 
the context of cuprates, one should however realize that 
this small- momentum peak structure in % q , hence also 
in j(k, q) is in fact suppressed by long-range Coulomb 
repulsion. Moreover, a bare matrix element for bond- 
stretching phonons is itself vanishing as q 2 at small q [ob- 
serve the matrix elements sin(g Q a/2) in Ea. (|43|) and re- 
sulting form-factor (1— 7 q ) in Ea.l|45|)]. These two factors 
eliminate the "forward-scattering" feature in cuprates. 
Altogether, it seems that effective coupling of the Fermi- 
surface electrons to bond-stretching phonons is somewhat 
reduced from A' -* for both small and large values of 
momentum-transfer q. Therefore, a significance of the 
bond-strcching phonons for the electronic properties of 
doped cuprates should not be overemphasized^ 



VII. CONCLUSIONS 

We calculated a momentum dependence of the static 
charge susceptibility in the t — t' — J model at various 
doping levels. We employed the formalism, which effec- 
tively resums the RPA series in Gutzwiller, J-term and 
Coulomb repulsion V^-channel simultaneously. We ob- 
serve that x(q) is a featureless function in the limit of 
weakly correlated overdoped regime. With decreasing 
of the doping level, strong correlations lead to a non- 
trivial, highly momentum-dependent renormalization ef- 
fects that cannot be described in terms of simple hole- 
density dilution effects. We demonstrate that x(q) is 
strongly suppressed near (tt,tt). However, x(q) remains 
large around (7r,0) and (0,7r) regions even at dopings as 
small as 0.1-0.2. A strong anisotropy of charge dynamics 
at finite wave-vectors has important experimental con- 
sequences, leading, e.g. to a pronounced anisotropy in 
a renormalization of the bond-stretching phonon modes 
as observed in cuprates. We find that the exchange 
J and second hopping t' influence the charge suscep- 



tibility in the opposite way - the former enhances it 
while t' and also Coulomb interactions weaken the small- 
momenta anomalies in x(q)- Implications of these find- 
ings on the spin and electronic properties deserve further 
analysis. The present study specifies the most "danger- 
ous" regions in a momentum space where one may expect 
charge-related anomalies. A complex interplay between 
the charge, spin and fermionic excitations should be con- 
sidered to fuller extent in order to locate more precisely 
a momentum position of low-energy charge modulations. 
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APPENDIX: COULOMB INTERACTION ON 
THE LATTICE: POINT-CHARGE LIMIT 

Consider the formula for the Coulomb interaction i|35|) 
for /(r) = 8(r). Introducing an auxiliary integration, the 
Fourier transform Vq at q z = can be written as 



Va 



V 

73/2 



POO 

J ° r^O 



2/2 , ~2 



(A.l) 



where rj 



X; 



yf denote square lattice sites with 
Xi, Vi = (0, ±1, ±2, . . .), while Zi = (0, ±1, ±2, . . ,)d ; we 
set a = 1 here. We use now the definition of the Jacobi 
i9 function, 



2 2iun 



n— — oo 



and represent Eq. I|A.1|) in the form 



V a 



V 



dr 



(A.2) 
(A.3) 



This formula can be considered as the limiting case of 
Eq. (|38|) with k — > oo. Numerically, the values of Vq in 
the whole Brillouin zone given by Eqs. and if^T^I 
become very close at n > 3. 
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